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The CHORUS experiment was designed to search for — > v r oscillation by detecting the 
decay topology of the r in charged current (CC) v T events. The detector was exposed to the 
Wide Band Neutrino Beam of the CERN SPS during the years 1994-1997. About 10% CC 
events were collected in the nuclear emulsion target. Up to now, about 170,000 events have 
been located and analysed. The speed of the automated emulsion scanning systems increases 
each year. With the present performance of these systems, it has become possible to perform 
large volume scanning. All tracks belonging to an interaction vertex can be recognized and 
measured precisely. This technique is not only applied to the search for neutrino oscillation 
but can also be used for the recognition of events where charmed particles are produced. 
Results obtained from the analysis of a sub-sample of the data on the production rate in 
CC interactions of neutral charmed mesons (D°) and charmed baryons (A c ) are presented. 
In addition a new measurement of the branching ratio for the decays of charmed hadrons 
into muons is given. Also measurements of topological branching ratios of D° and A c are 
presented. Finally, a search for associated charm production is discussed. 



1 Introduction 

1.1 The CHORUS experiment 

The CHORUS experiment was proposed primarily to search for Vn — > v T oscillations through 
the appearance of v T in a beam, aiming to explore the domain of small mixing angles down 
to sin20^ r ~3x 10~ 4 for mass parameters Am 2 ~ 1 eV 2 . This represents an order of magnitude 
improvement over the previous generation of experiment^. 

The CHORUS detector is a hybrid set-up that combines nuclear emulsion stacks with various 
electronic detectors. The nuclear emulsion acts both as target for neutrino interactions and as 
detector, allowing three-dimensional reconstruction of short-lived particles as the r~ lepton and 
charmed hadrons. The nuclear emulsion target, which is segmented into four stacks, has an 



overall mass of 770 kg, each of the stacks consisting of eight modules of 36 plates of size 36 x 72 
cm 2 . Each plate has a 90 /im plastic support coated on both sides with 350 /im emulsion layers. 
Each stack is followed by three interface emulsion sheets having a 90 /im emulsion layer on 
both sides of an 800 /mi thick plastic base and by a set of scintillating fiber tracker planes. 
The interface sheets and the fiber trackers provide accurate particle trajectory predictions into 
the emulsion stack in order to locate the vertex positions. The accuracy of the fiber tracker 
prediction is about 150 /im in position and 2 mrad in the track angle. The major drawback is 
the absence of the time information: any charged particle traversing the emulsion during the 
exposure will leave a track. 

The emulsion scanning is performed by fully automatic microscopes equipped with CCD 
cameras and a read out system, called Track Selector*. In order to recognize track segments 
in an emulsion, a series of tomographic images are taken by focusing at different depths in the 
emulsion thickness. The digitized images are shifted according to the predicted track angle and 
then added. The presence of aligned grains forming a track is detected as a local peak of the 
gray level of the summed image. The track finding efficiency of the track selector is higher than 
98% for track slopes less than 400 mrad. 

The electronic detectors downstream of the emulsion target include a hadron spectrometer 
which measures the bending of charged particles in an air-core magnet, a calorimeter where the 
energy and direction of showers are measured and a muon spectrometer which measures the 
charge and momentum of muons. 

1.2 Data Collection 

The West Area Neutrino Facility (WANF) at CERN provides a beam of 27 GeV average energy 
consisting mainly of with a 5% contamination. For the four years in which the emulsion 
target was exposed, the integrated beam intensity corresponds to 5.06xl0 19 protons on target. 
The analysis of the data from the electronic detectors allows the identification of the set of events 
possibly orginating from the emulsion stacks. For the first phase of the analysis, the events were 
subdivided into two classes, based on electronic detectors, which are the one-muon and the zero- 
muon samples, distinguished by the presence or absence of one reconstructed muon of negative 
charge. For vertex location, track trajectories belonging to reconstructed events in the electronic 
detectors are used to guide the scanning. The track is first searched in the interface sheets and 
then followed back into the target emulsion. If it is not found in two consecutive plates, the 
first of these is defined as vertex 'plate. This plate may contain the primary vertex or the decay 
vertex, or both. Once the vertex is located, the charm decay search is performed using the 
netscan method^. It consists of recording all track segments within an angular acceptance in a 
volume surrounding the assumed vertex position. 



2 Charm analysis 

Charm production in neutrino charged current (CC) interactions has been studied in several ex- 
periments^ with electronic detectors and mostly through the analysis of dimuon events. In these 
events, the leading muon is interpreted as originating from the neutrino vertex and the other, of 
opposite charge, as the product of the charmed particle semileptonic decay. Experiments of this 
type, however, suffer from a significant background of events (~ 30%) in which the second muon 
originates from an undetected decay in flight of a pion or a kaon rather than from a charm decay. 
Moreover, the type of the charmed particle and its decay topology can not be identified in these 
experiments. Nevertheless they have provided measurements of the strange quark content of the 
nucleon as well as an estimation of the charm quark mass. A much lower level of background 
can be achieved using an emulsion target which provides a sub-micron spatial resolution, and 



hence, the topological identification of charmed hadron decays. The statistics accumulated in 
this way is however very limited. Only recently, the development within CHORUS of automatic 
scanning devices of much higher speed has made studies of charm production with high statistics 
possible. 

2.1 D° production rate measurement 

In the past, D° production in neutrino interactions was studied in E531 experiment. Only 
57 D° events were identified and analysed. In CHORUS, we have performed an analysis with 
much higher statistics. About 25,000 one-muon events were analysed representing ~15% of the 
CHORUS data. The main criteria for the event selection are the following: at least one of the 
decay daughters and muon track reconstructed with more than one-segment, matched tracks in 
the electronic detector; and the daughter track must have a significant impact parameter with 
respect to the vertex point. These criteria select 851 events from 25,000 CC Vn events. The 
selected events are visually inspected in order to distinguish decays from hadronic interactions, 
gamma conversion and false vertices which are reconstructed using background tracks. After 
performing visual inspection, we confirmed that out of 851 selected events, 226 show 2-prong 
decay topology and 57 show 4-prong decay topology. 

The efficiency of the charm selection was evaluated with a GEANT3 based simulation of 
the experiment. Large samples of deep-inelastic neutrino interactions were generated according 
to the beam spectrum using the JETTA generator derived from LEPTO and JETSET. The 
simulated response of the CHORUS electronic detectors is processed through the same recon- 
struction program used for data. To evaluate the netscan efficiency, realistic conditions of track 
densities need to be reproduced. This is achieved by merging the emulsion data of the simulated 
events with real netscan data which do not have a reconstructed vertex but contain tracks which 
stop or pass through fiducial volume representing the real backgrouncP. The combined data are 
passed through the same netscan reconstruction and selection programs as used for data. The 
selection efficiencies are estimated as (58.6±5.1)% and (70.1±5.2)% for 2-prong and 4-prong de- 
cays respectively. Based upon 282 D° decays with an estimated background of 9.2±2.6 K° and 
A , we obtaiiP ff(J " g "ff V4) =(1.99±0.13(afafc)±0.14 (syst.)) xl0~ 2 at 27 GeV. The topological 
ratio V4/V2 is found to be 

(23.1 ±4.0)% (1) 

Combining (1) with Br(D° -» U4) = (13.3 ± 07) x 10~ 2 using the PDG tabled, we can obtain 
branching ratio of D decaying into neutral particles: 

Br (D° - neutral) = 1 - Br(D° U4) x (1 + ( |^g^Z| ( 2 ) 



= (29.1 ± 10.4) x 10~ 2 
The precision of this measurement will be improved with the final statistics. 



2.2 Semi-leptonic branching rate measurement 

To measure the semi-leptonic branching ratio, the event the event selection criteria are the 
following: at least two tracks of more than one-segment coming from different vertices match 
tracks at the electronic detectors. About 50,000 CC events are analysed with this selection 
and in total, 1055 events are selected. To estimate the selection purity, 25% of selected events 
are visually inspected. The selection purity is obtained as 91%. Based on this selection purity, 
the number of charmed hadrons is 956±35. 



The efficiency of charm selection was evaluated as explained in the previous section. The 
correction factor (R) for the efficiency in the determination of the semileptonic branching fraction 
is defined as where e D i is the selection efficiency for charm species Di , e D v the selection 

i ^ i 

efficiency for semileptonic decays of Di, the fragmentation fraction. Based on the estimated 
efficiencies R is l.Olrfc 0.05. 

The average semi-leptonic branching fraction can be written in terms of measurable quanti- 

ties as = N 2 £ L x R where N seL is the number of selected events, corrected for the selection 
purity and N^' is the number of selected events with a secondary muon in the final state 
corrected for selection purity as well as the muon identification efficiency and purity. 
Combining these numbers, we measured = 0.0093 ± (stat.) ± 0.009{syst.j& 

2.3 A c production rate measurement 

Although evidence for charmed baryon production by neutrinos has been reported in the litera- 
ture with few events observed in a number of different bubble chamber experiments, cross-section 
values are known with a large error. 

Our analysis is based on a statistical approach using the flight length distributions of charmed 
hadrons. Since the lifetime of the A c is smaller than that of other charged charm hadrons, the 
sample of charged charm decaying at distances less than 400 fj,m from the v vertex, called short 
decays, should be dominated by A c . Conversely, long decays are enriched by D + and Df decays. 
Therefore, the analysis is performed applying two different selections. Event selection for short 
flight decays is based on the following: muon track and daughter track must be reconstructed 
having more than one segment and the daughter track must have a big impact parameter with 
respect to the vertex point (5/xm < IP < 30 /mi). We have analysed about 50,000 CC v events 
with short decay selection criteria. In total, 1614 events were selected for visual inspection. 62 
events are confirmed as 1-prong decay and 66 events show 3-prong decay topology. 

On the other hand, for long flight decays we require that the parent track and muon track 
must be reconstructed by the reconstruction software and the parent track must have a small 
impact parameter with respect to the neutrino vertex and at least one of the daughter tracks 
must be reconstructed with more than one segment and minimum distance between parent 
track and daughter track must be less than 5 /xm. About 56,000 CC v events were analysed and 
586 events were visually inspected. The confirmed charm events after flight length cut are the 
following: 133 1-prong and 195 3-prong decays. 

Combining short and long decay search and taking into account efficiency and background, 
the number of A c is 861±198(stat.)±98(sysf.)jl 5 4 (QE). Based on these events the branching 
ratio into 3-prongs was determined: BR(A C — ► 3 pron<7)=0.24±0.07(stat)±0.04(syst.) and the 
total charged-current cross sectioM f^^=(1.54±0.35(stoi.)±0.18 (syst.)) x 10~ 2 . 

2.4 Associated charm production 

Associated charm (cc) production in neutrino interactions is a very rare process and therefore 
very difficult to observe. In the past, cc production in charged-current Vn interactions has 
been studied through trimuons and like-sign dimuons , however, the reported production rate is 
higher than theoretical expectation. On the hand, only one event consistent with neutral-current 
(NC) production of cc has been observed by the E531 collaboratioffl Based on one event, the 
production rate with respect to NC production is estimated as =(0.13±0.31) xl0~ 2 . The 
sub-micron resolution of nuclear emulsion allows the study of this kind of rare processes. We 
have performed a search in CC interactions and one event is observed with the characteristics 
of cc production?^. A new search has been started, both in the NC and CC event sample. Fig. 
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Figure 1: Sketch of the event topology in emulsion. 

1. shows the schematic view of one of the candidate events. This event is found in the zero-muon 
sample and shows two V2 decay topologies. One minimum-ionizing particle and three heavy 
ionizing particles are emitted from primary vertex. A neutral particle(C°) decays in the same 
emulsion plate, 62 /iin downstream of the primary vertex. In the next plate, another neutral 
decay is found with the flight length of 9 77 fim. The lifetimes of the decaying particles are 
estimated using the relation: T= l<i>\n\ where < 9 > is average emission angle of daughter 
particles relative to the parent direction , I is the flight lenght of the decaying particle and c 
is the speed of light. The estimated lifetimes, shown in table 1, are consistent with the decay 
of the neutral charmed particles. The search for new candidates and more detailed analysis of 
already observed events are in progress. 



Table 1: List of particles at primary and secondary vertices. 





9 y {rad) 


6 z (rad) 


l(pm) 


r(s) 


particlel 


-0.060 


-0.002 






1-C U 


0.008 


0.031 


63 


2.6 x 10" 14 


daughter (at 1-C U ) 


-0.019 


0.040 






daughter (at 1-C U ) 


-0.052 


0.028 






2-C u 


-0.078 


0.095 


977 


5.5 x lO" 13 


daughter (at 2-C u ) 


-0.166 


0.049 






daughter (at 2-C u ) 


0.035 


0.103 
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